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ABSTRACT 



This is the third paper in a series in which we attempt to put constraints on the flattening of dark halos in disk galaxies. We observed 
for this purpose the HI in edge-on galaxies, where it is in principle possible to measure the force field in the halo vertically and radially 
from gas layer flaring and rotation curve decomposition respectively. For this purpose we need to analyse the observed XV diagrams 
in such a way as to accurately measure all three functions that describe the planar kinematics and distribution of a galaxy: the radial 
HI surface density, the rotation curve and the HI velocity dispersion. In this paper, we first present the results of the modelling of 
our HI observations of 8 HI rich, late-type, edge-on galaxies. We find that in all of these we find differential rotation. Most systems 
display HI velocity dispersions of 6.5 to 7.5 km s -1 and all except one show radial structure in this property. There is an increase in 
, the mean HI velocity dispersion with maximum rotation velocity, at least up to 120 km s _1 . 

Next we analyse the HI observations to derive the radial variation of the thickness of the HI layer. The combination of these gas flaring 
measurements with the HI kinematics measurements allow us to measure the total vertical force of each galaxy assuming hydrostatic 
equilibrium. We find that with the exception of the asymmetric IC5052, all of the galaxies in our sample are good candidates for 3D 
mass modelling to measure the dark halo shape. The flaring profiles are symmetric with respect to the galactic centres and have a 
common shape, increasing linearly inside the stellar disks and exponential outside where the gravitational potential is dominated by 
the dark halo. In the best example, UGC7321, we find in the inner regions small deviations from the midplane and accompanying 
increases in thickness of the HI layer that are possibly a result of perturbations of the gravitational field by a relatively strong bar. 
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> ' 

ON 1 . Introduction force field in the vertical direction for a Gaussian gas layer can 

be written as 

In paper I in this series (O'Bri en et alll20 10a) we presented HI 



', observations of a sample of 8 edge-on, HI rich, late-type galax- dK z c"v,g 



(1) 



£f) ■ ies. The aim of the project has been described there in detail. Q z ' (FWHIVL ? /2.35) 2 ' 
<^> Briefly, we attempt to put constraints on the flattening of dark ha- 

<^> , los around disk galaxies by measuring the force field of the halo where cr 2 g is the vertical velocity dispersion of the gas and 

,— 1 1 vertically from the flaring of the HI layer and radially from ro- FWHIVL ^ the gas layer thickness. Using the flaring of the HI 

\ tation curve decomposition. For the vertical force field we need distribution and the derived density distribution of the gas and 

• i-H i to determine in these galaxies both the velocity dispersion of the stellar distributions, it is possible to measure the halo shape over 

^ ' HI gas (preferably as a function of height from the central plane the HI extent of the luminous disk. For a given vertical gas ve- 

H ] of the disk) and the thickness of the HI layer, all of these as a locity dispersion, the flaring will decrease and the gas surface 

. . . function of galactocentric radius. In addition we also need to ex- density will increase, for increasingly flattened dark halos. Thus, 

tract information on the rotation curve of the galaxy and the de- in gas-rich late-type galaxies it is possible to measure the force 

projected HI surface density, also as a function of galactocentric field of the halo vertically from the gas layer flaring and radially 

radius. To adequately constrain the halo shape these parameters from rotation curve dec omposition. 

needed to be measured out to as low surface densities as possi- In the second paper (lO'Brien et all 1201 Obi, paper II, we pre- 

ble.Here we present the determination of these properties from sented a new method to accurately determine the rotation curve, 

our HI observations. deprojected HI surface density and velocity dispersion at all radii 

In order to resolve the vertical structure, we selected a sam- in a edge-on gas disk. The superposition of velocity profiles from 

pie of nearby, HI rich, late-type edge-on galaxies which are man y radii in each sightline through an edge-on HI disk tends to 

likely to exhibit large HI flaring expected from high dark mat- cause an overestimate of the velocity dispersion with most mea- 

ter fractions. From paper I, we recall that the gradient of the surement methods (see paper II). Measuring the radial flaring 

profile requires a model of the galaxy rotation and face-on sur- 
face density; this also necessitates high accuracy rotation curve 



For correspondence contact Ken Freeman or Piet van der Kruit measurement and HI surface density deprojection. 
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Fig. 1. ESO074-G015 (IC5052). This figure shows the results of the modelling of the kinematics. The best fit measurements for 
each kinematic function are shown in as (red) crosses for the the HI velocity dispersion in the top left-hand panel and in the lower 
left-hand plot for the deprojected HI surface density. The top right-hand panel shows the derived rotation curve (dark/red) overlaid 
on the observed XV diagram. Also shown on the XV diagram are the rotation fits obtained using the less accurate envelope tracing 
method (dark grey/orange) and the peak flux method (light grey/yellow). The latter always shows a lower rotation velocity than the 
other two fits. The lower right-hand plot displays the derived rotation curve with the receding and approaching sides separately (in 
the electronic version red and blue respectively). The FWHM synthesised beam has dimensions 292 x 292 pc. 



In this third paper we present the rotation curve, HI veloc- 
ity dispersion and deprojected HI surface density from the XV 
maps measured using the iterative XV diagram modelling pro- 
gram (see paper II). Previously, measurement of these parame- 
ters in edge-on galaxies has been difficult. Indeed, the gas ve- 
locity dispersion has been measured as a function of radius in 
only one edge-on galaxy NGC4244 ( Ollingl IT996ah . If the data 
have suffiently high S/N, then our iterative XV modelling pro- 
gram is capable of deriving the variation of rotation, velocity 
dispersion and density as a function of both radius and height 
z above the equatorial plane of these edge-on galaxies. For the 
data described in paper I, the S/N is not high enough to derive 
the variations in both z and R. We use the XV diagrams inte- 



grated over z, and measure the gas surface density distribution 
and kinematics as a function of radius only, assuming that the 
gas velocity dispersion is isotropic and vertically isothermal. 



2. Results for kinematic modelling 

For each galaxy in our sample we present the best fit measure- 
ments of the radial kinematics, and show the map of the residuals 
formed by comparing the model XV diagram to the observed XV 
map. These XV maps were formed by integrating the observed 
HI channel maps over the z axis. 



J.C. O'Brien et al.: The dark matter halo shape of edge-on disk galaxies III 



3 




Radial surface profile Rotation curve 




-20 -10 10 20 5 10 15 

Radius (kpc) Radius (kpc) 



Fig. 2. ESO109-G021 (IC5249). This figure shows the results of the modelling of the kinematics. For an explanation of the figure 
see the caption of Fig.Q] The FWHM synthesised beam has dimensions 1 179 x 1 179 pc. 



2.1. ESO074-G015(IC5052) 

Although ESO074-G015 is one of the nearest edge-on low mass 
galaxies, it is a galaxy less suitable for dark halo fitting be- 
cause it is asymmetric in both the HI and optical distribution. 
Furthermore, it shows a bright star formation region 2 kpc west 
of the galactic center with associated radio continuum emission, 
indicating a possible regional deviation from hydrostatic equi- 
librium. The measured rotation curve (see Fig. [TJ also shows a 
large asymmetry which is unlikely to be solely due to the HI 
warp. Given the irregular optical distribution, it is more likely to 
indicate a non-axisymmetric mass distribution commonly found 
in Magellanic dwarf galaxies. The centre of mass was poorly 
determined in IC5052 due to the irregular structure of the XV 
distribution. 

Despite the apparent mass asymmetry, the HI velocity dis- 
persion is roughly constant at 6-7 km s _1 . Similarly the HI ra- 
dial profile is reasonably symmetric, although it displays a broad 



plateau in the bright star formation region on the west side. As 
the channel maps clearly show vertical flaring of the HI layer, we 
will proceed with the flaring analysis for this galaxy. However, 
the gross kinematic asymmetry precludes it from use in the hy- 
drostatic analysis to determine the shape of the dark matter dis- 
tribution. 

2.2. ESO109-G021 (IC5249) 

ESO109-G021 (IC5249) is one of several well-known "su- 
perthin" galaxies, renowned for the large optical axial ratio. The 
rotation curve shown in Fig. [2] shows that ESO109-G021 has 
a maximum rotation speed of 1 10 km s _1 . ESO109-G021 is the 
only one of the galaxies in our sample with rotation speeds larger 
than 100 km s" 1 that does not show the "figure-8" signature in 
the HI distribution indicative of a possible bar. 

The measured HI radial surface density has a shallow 
gradient with dips at around 11-12 kpc in agreement with 
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the features in th e depro jected surface density obtained by 
Ivan der Kruit et al. (2001). However these authors also mea- 
sured a central HI depression, while we find a small central peak 
when measuring from the left-hand side and a plateaued central 
surface density when measuring from the right-hand side. The 
deprojected surface density was measured after folding the ob- 
served surface density calculated from the XV map about the 
galactic center. 

One possible explanation for th e discrepancy between ou r 
measured surface density and that of Ivan der Kruit et"aT] (1200 lh . 
is the lower sensitivity of their observations. In addition to the 
12 hr 6.0C ATCA observation used by all parties, we acquired 
an additional 30 hr of HI synthesis observations spread over the 
ATCA EW352, 750D and 6.0A array configurations. This re- 
sulted in a peak signal-to-noise of 46.1 and a synthesised beam 
FWHMe of 8" in the XV map, which appears to be a signifi- 
ca ntly higher sensitivity and spatial resolution than that acquired 
bv lvan der Kruit et al.1 il200il) . 



The measured rotation curve is reasonab ly symmetrical , 
in agreement with that obtained by Ivan der Kruit et al .1 (2001) 
which was derived by modelling the XV distribution us- 
ing the fixed veloc ity dispersion XV modelling method of 
van der Kruit ( 1981, see also Sect. 2. 5 in paper II). Like that of 
van der K ruit et al.l(l200ll) . our measurement shows a steep inner 
rise — as oppos ed to the earlier less accurate measurements by 
lAbe et al.l d 19991) which appear to have been measured from the 
peak flux ridge of the XV map. 



The measured HI velocity dispersion varies from 5.5 to 8.5 
km s _I , falling to lower values on the right-hand side of the 
galaxy. The mean HI velocity dispersion of 6.9 + 0.8 km s is 
in agreement with the b est fit value of 7 km s _I found by 
Ivan der Kruit et all (l200lh . 
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Fig. 4. ESO138-G014. This figure shows the results of the modelling of the kinematics. For an explanation of the figure see the 
caption of Fig.Q] The FWHM synthesised beam has dimensions 965 x 965 pc. 



2.3. ESO115-G021 

The high resolution HI imaging of ESOH5-G21 yields an al- 
most solid-body rotation curve flattening out at a radius of 5-6 
kpc at a peak rotation speed of approximately 65 km s _1 (see 
Fig. [3). Like the other small galaxies in our sample it shows 
steep inner rotation unlike the solid body rotation commonly ob- 
tained from low resolution studies of small galaxies. Despite an 
apparently accurate rotation curve fit there is still significant flux 
in the residual in the part of the XV diagram containing flux at 
large radii. This appears to be due to an underestimate of the HI 
surface density and the velocity dispersion at large radii. 

The measured HI velocity dispersion is surprisingly low at 
only 4.5-6 km s _1 at radii within 5 kpc from the galactic cen- 
tre. This galaxy displays particularity bright HI emission in the 
inner 2-3 kpc, with a brightness temperature Tb ~ 150 K. The 
low measured velocity dispersion and high brightness tempera- 
ture are consistent with a partially opaque inner HI disk. In the 



face-on galaxy NGC1058 jPetric & Rupeiil (f2Q06) found regions 
of HI velocity dispersion as low as 4 km s _I — indicating cold 
low velocity dispersion HI on scales of at least 500 pc which 
is quite large compared to typical star formation scales in the 
ISM. Given the large beamwidth of 1.5kpc of the NGC1058 HI 
observations, it is likely that smaller regions of even lower veloc- 
ity dispersion gas exist in NGC1058. In ESOl 15-G021, regions 
with a HI velocity dispersion of 4.5 km s _1 occur on radial scales 
of 1 kpc, similar to the cold gas distribution in NGC1058. 

2.4. ESO138-G014 

ESO138-G014 is a relatively large Sd galaxy with a HI disk ex- 
tending to 20 kpc, and a maximum rotation speed of 112 km 
s (see Fig. @J. The rotation curve is clearly differential, rising 
steeply in the innermost 5 kpc, then starts to flatten outside of 10 
kpc from the galactic centre. As noted in paper I, the XV map 
displays a faint "figure-8" signature, indicating a possible bar. 
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Fig. 5. ESO146-G014. This figure shows the results of the modelling of the kinematics. For an explanation of the figure see the 
caption of Fig.Q] The FWHM synthesised beam has dimensions 793 x 793 pc. 



The HI velocity dispersion is approximately 7 km s at radii 
outside of 9 kpc. At radii inside of 9 kpc, the HI velocity disper- 
sion rises steeply to a plateau of 10-11 km s _1 within 5 kpc radii. 

The HI radial profile is roughly constant within ~9 kpc from 
the centre, except for central peak within ~2 kpc radius. The 
large central plateau spans the same radial domain as the in- 
creased velocity dispersion and could be explained by radial 
streaming in the bar causing redistribution of the gas. As ex- 
plained in Sect. 14.41 we will leave this galaxy out of the discus- 
sion because of missing short spacings and a poorly determined 
flaring profile of the HI layer. 

2.5. ESO146-G014 

ESO146-G014 is another small Sd galaxy similar to ES0115- 
G021 with a particularly low HI velocity dispersion of ~4.5 km 
s _1 at radii less than 4 kpc. At larger radii, the mean velocity 
dispersion is 6 km s _1 , with a peak of 7.5 km s _1 at 5 kpc west 



of the galactic centre. This high velocity dispersion region on 
the west side corresponds to a 1 kpc wide bright overdensity in 
the B-band optical image, and to a 1 kpc gap between bright HI 
regions in the HI total intensity map (see paper I, Fig. 5a). 

The rotation curve of ESO146-G014 peaks at only 70 km 
s _1 and is quite symmetric about the galactic centre, particularly 
within radii of 4.5 kpc. The HI surface density shows a bright 
central peak, a plateau from 3 to 6 kpc radius, after which it 
declines steeply on the left-hand side, and less so on the right- 
hand side. 

2.6. ESO274-G001 

ESO274-G001 is the closest galaxy in our sample at a distance 
of only 3.4 Mpc. This allowed the 6 km maximum baseline of 
the ATCA to yield an exceptionally small HI synthesised beam 
of 162 pc. Thus despite the small radius of only 7-8 kpc, each 
side of the galaxy is spanned by over 47 independent beams. 
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Using a beam sampling of 3 pixels per beam, approximately 150 
kinematic fits are obtained on each side of the galactic centre. 



The measured kinematics shown in Fig. [6] display a differ- 
entially rotating galaxy with a maximum rotation speed of only 
86 km s _1 . This rotation velocity shows that ESO274-G001 is a 
small galaxy with a rotation curve that continues to rise until the 
outermost kpc. The 20%-level envelope fitting method overesti- 
mates the rotation speed by up to 10-15 km s _1 , while the peak 
flux method dramatically underestimates it by up to 50%. The 
HI surface density is reasonably symmetric, with a central peak 
within 2 kpc, and a plateau to 5.5 kpc on the left-hand side and 
4.5 kpc on the right-hand side, declining to an HI surface density 
of 7-8xl0 19 atoms cm . The HI velocity dispersion at 6.5 km 
s _1 shows no systematic trend with radius, with small oscillations 
of only 1 km s _1 around this value, throughout the HI disk. 



2.7. ESO435-G025 (IC2531) 

ESO435-G025 is the only Sc galaxy in our sample and more 
massive than the next largest galaxy by a factor of four. The 
rotation curve displays an unusual concave shape (see Fig. [7]) 
indicating a dip in the gradient of rotation. ESO435-G025 
was one of the first dete ctions of a bar in an edge-on galaxy 
dBureau & Freemanlll999l) . Despite the bar, the measured kine- 
matics are an excellent representation of the galaxy. The inner 
rotation of IC253 1 is similar to that of NGC891, however lack- 
ing g as tracing an inner ring, seen in both HI dSwaters et al.L 
1997) and CO (Tony Wong, private communication). It is 
quite possible that like NGC891, the inner galaxy of ES0435- 
G025 is largely populated by mole cular gas in the inner disk 
dWong. Blitz & van der HuTsll2004l) . 

ESO435-G025 is the most distant galaxy in our sample and 
with a low peak signal-to-noise of 41.4 in the XV map, and par- 
ticularly low signal-to-noise in the parts of the XV diagram with 
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Fig. 7. ESO435-G025 (IC2531). This figure shows the results of the modelling of the kinematics. For an explanation of the figure 
see the caption of Fig.Q] The FWHM synthesised beam has dimensions 1305 x 1305 pc. 



radii less than 7 kpc. Unfortunately, the signal-to-noise in the 
inner galaxy was too low to deconvolve for the effect of the HI 
beam (FWHMg = 1.3 kpc) while modelling the XV diagram. 
As a result the HI surface density is probably overestimated at 
radii within 10 kpc where the gradient of rotation over each 
beamwidth is large. The measured HI velocity dispersion dis- 
plays a significant peak and central trough, which is largely an 
effect of beam smearing. At radii outside 10 kpc, where the beam 
smearing effect is minimal, the HI velocity dispersion has a mean 
of 7 km s" 1 . To compensate for the inability to properly measure 
the velocity dispersion inside of 10 kpc, we adopt a smooth cen- 
tral velocity dispersion rising from 8 km s _1 at 10 kpc radius 
to 10 km s _1 in the centre. Inside of 10 kpc we adopt a central 
plateau for the HI surface density, but we flag that the derived 
properties within 10 kpc should be treated with caution. 



2.8. UGC7321 

The HI bar of UGC7321 causes a significant distortion in the 
peak flux of the XV distribution. The effect of non-circular mo- 
tions in the bar causes a deviation in the projected rotation of gas 
rotating in the bar. The XV modelling is able to recover the az- 
imuthally averaged surface density, and rotation (Fig. [8a]>; how- 
ever, it is unable to fit the non-axisymmetric gas in the bar, as 
can been appreciated from Fig. [8b] which shows the observed, 
model and residual XV diagrams. Orbits in a bar r ed po tential 
have been studied by IContopoulos & Mertzanidei d!977l) . who 
found that there are two families, where the main orbit family 
is designated x\ (see also lBinnev et all 1 199 lh . On the left-hand 
side of the galactic centre, the X\ orbit showing the near and far 
side of the bar is clearly apparent in the residual map. The resid- 
ual structure on the right-hand side of the galactic centre is less 
clear. 

The peak flux in the x\ orbits of the residual is a 15% of 
the peak flux of the observed XV diagram. To ascertain the 
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Fig. 8a. UGC7321. This figure shows the results of the modelling of the kinematics. For an explanation of the figure see the caption 
of Fig.Q] The FWHM synthesised beam has dimensions 727 x 727 pc. 



likely surface mass density perturbation caused by the bar, we 
inspected the K band azimuthal flux profile of the bars in several 
late-type face-on galaxies with weak bars. If one assumes a sig- 
nificant fraction of the brightness variation is biased by young 
stars contributing to the bar, galaxies with a relatively weak, dif- 
fuse bar are unlikely to have a difference of stellar surface den- 
sity between bar and inter-bar regions greater than 50% (for ex- 
ample, the Sc galax y NGC2525). Photometry of UGC7321 by 
iPohlen et al.l {2003) shows that the stellar bar is very faint and 
thin, indicating that minimal buckling of the bar has occured. 
The thin stellar bar of UGC7321, and the small HI mass frac- 
tion in the bar suggest that the inner mass density perturbation 
is quite low, and unlikely to significantly effect the halo mass 
modelling. 



3. The flaring of the Hi layer 

There have been few attempts to directly measure the flaring of 
the HI laye r in galaxies. W ith the exception of Oiling 's study of 
NGC4244 (Oiling, 1996ah and a n earlier less sensitive study of 
NGC5023 (Bottem a etail 1 19861) . all the flaring measurements 
were undertaken on massive Sb-Sc galaxies. It is particularly 
difficult to measure the flaring in large galaxies, as such galax- 
ies are expected to have higher disk surface densities and conse- 
quently thinner gas layers. For gas in hydrostatic equilibrium the 
FWHM,(/?) of the HI is inversely pro portional to th e maximum 
rotation velocity V max (Eqn. (D31) of foiling! [l995h 

FWHM z hi{R) oc y A \ (2) 

V max 

so large galaxies with max rotation speeds similar to the Milky 
Way, will only flare to 1/3 of the height of small galaxies like 
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Fig. 8b. UGC7321. This figure shows the observed XV diagram (upper panel), the model XV diagram formed from the best fit 
kinematics (middle panel), and the map with residuals between the two (lower panel). The contours are in units of the rms noise, 
which is <t = 0.524 mJy beam -1 . The FWHM of the HI synthesised beam is 727 pc. To the left of the galactic centre the residual 
map clearly shows the gas traced by the orbits in the bar. This is less clearly seen on the right-hand side due to the small warp. If 
one excludes the gas in the bar, the fitted kinematics depict a very good representation of the observed XV map. 



ESO146-G014, given similar HI vertical velocity dispersions 
o~v,hi(R)- Most of these previous HI flaring measurements of 
large galaxies were undertaken around 20 years ago, with larger 
telescope beams and poorer sensitivity resulting in partial reso- 
lution of the FWHM, vertical thickness of the gas layer. 

One of the first measurements of the gas layer thickness in 
an external galaxy was the early lo w resolution WSRT study of 
N GC891 bvlSancisi & Aiienid 19791) . This data were re-analyzed 
by Ivan der Kruitl (1198 lh - These authors found that the vertical 
thickness of the HI layer was unresolved by the 1 .4 kpc beam 
in the inner parts (R < 8 kpc) and increased to 2 kpc at the 
outermost measured point of R » 14 kpc. Using deeper higher 
resolution VLA data with synthesised beam FWHMe of 0.4 and 



0.9 kpc respectively for NGC4565 and NGC891. lRuped dl989h 
found that the gas thickness increased exponentially with radius, 
although NGC891 displays substantial scatter. This is due to a 
low channel map peak signal-to-noise of « 13 in the NGC891 
HI cube, compared to as 64 in the NGC4565 dataset. This result 
is consistent with the thickness of the gas layer in our Galaxy, 
which is also found to increase exponentiall y from 200 pc nea r 
the Galacti c Centre to 1 .3 kpc at -20 kpc (iMerrifieldl 1 19921) . 
Meanwhile, Brau fitted the gas thickness of M31 with an 

exponential z profile and found that the gas thickness increases 
roughly linearly from a scale height of 200 pc at 5 kpc radius to 
600 pc at 25 kpc. 
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The most accurate measurement of HI flaring w as obtained 
for the small Scd galaxy NGC4244 by Oiling ( 19963) who found 
that the FWHM, thickness increased from 400 pc at 5 kpc radius 
to 1.8 kpc at 14 kpc radius, with uncertainty ranging from « 100 
pc at inner radii to 200-500 pc at the outermost measured point. 
These were by far the most sensitive high resolution observa- 
tions obtained, with 14 hrs of VLA data spread over B, C and D 
arrays, resulting in a FWHMg of 170 pc. The only other small 
galaxy that has been observed is the Scd NGC 5023 imaged at 
the W SRT with a spatial resolution of ~650 pc (Bott ema et all 
fT986h . These authors find a constant thickness of 1.2 + 0.3 kpc 
despite the fact that their Fig. 4 appears to show an increase 
of thickness from 1.0 to 1.7 kpc, particularly on the east side. 
However the spectral resolution of 16.6 km s _1 available on the 
WSRT at the time was a little coarse given the low maximum 
rotation speed (86 km s -1 ) of NGC5023. This low spectral res- 
olution and the low signal-to-noise of the observations would 
have increased the uncertainty in their radial flaring profile. 



thickness and flaring. This study also involved an investigation 
into the effects of small residual inclinations away from purely 
edge-on. 

The technique used bv lBraunl ([T99 1 ) to measure the flaring of 
the inclined galaxy M3 1 was more complicated. Braun used the 
fitted inclination and rotation curve to fit the vertical thickness 
of HI features in the spiral arms. 

In the next section, we present the method which we have 
developed to accurately deproject the projected radius to the 
galactocentric radius for each vertical thickness measurement, 
and correct the measured thickess for the spatial resolution of 
the telescope along the z axis. As our method fits the vertical gas 
thickness in all channel maps with sufficient flux, the number 
of measurements at each radius produces a more reliable radial 
flaring profile, and also allows the flaring profile to be measured 
at radii much closer to the ga l actic ce nter. Our method is very 
similar to that used by Oiling (1996a); however, we only allow 
for warps of the HI disk in the plane of the sky. 



3.1. Methods used to measure the vertical gas thickness 

Several methods have been used to measure the vertical flaring 
of high inclination galaxies. The gas disk thickness must be mea- 
sured from individual channel maps, rather than from the total HI 
column density map, as the latter measures the thickness of the 
integrated flux. If the HI distribution is well resolved by the tele- 
scope beam, measuring the gas thickness from channel maps is 
elementary. The challenge lies in deprojecting the measurements 
from the projected radius of each measurement to the galac- 
tocentric radius. In each method below, the measured FWHM, 
vertical thicknesses are corrected for the spatial resolution of the 
beam (see Eqn. (O below). 

The most common method used to determine the radial flar- 
ing profile is to fit the vertical thickness of just the extreme ve- 
locity channel maps along each line of sight. This provides an 
accurate measurement of the flaring in the outer disk if the ro- 
tation curve is flat over this radial domain. If gas in the extreme 
velocity channel maps is emitted from a narrow azimuthal range 
around the line of nodes, then the gas fitted by a vertical profile 
at each projected radius R' is emitted from a narrow radial range, 
allowing the radius of each thickness measurement to be ap- 
proximated by the major axis distance R'. This m ethod was used 
in the early fl a ring m easurements of NGC891 dSancisi & Alle"nl 
I1979L lRur)enLl l989h NGC4565 dRunenL Il989h and NGC5023 
dBottema et allll986l) . 

Unfortunately this deprojection approximation can often 
miscalculate the radial shape of the flaring profile. This occurs 
in less differentially rotating galaxies, as the radial range of gas 
emission is larger along each line of sight through the extreme 
channel maps. Also, less massive galaxies with shallower rota- 
tion curves will have a smaller major axis extent spanned in the 
outer velocity channels and the radial domain where the rotation 
curve is most flattened may not be where the gas disk flares. The 
other disadvantage is that the vertical thickness at each radius is 
measured only twice, once on either side of the galaxy. Another 
problem is that the gas thickness is only measured at radii with 
flux in the edge channel maps. In a large galaxy with a flat ro- 
tation curve, these channel maps will probe most of the radial 
extent of the galaxy, however in a less massive galaxy with a 
less steep inner rotation curve, the radial span of the outer chan- 
nel maps is much less. 

The analysis oflvan der Kruitl (Il98ll) of the NGC 89 1 data of 
ISancisi & Alien! ( 1 19791) involves matching the observed XV map 
to a simulated XV formed with different assumed values of gas 



3.2. Gas flaring measurement by iterative HI cube modelling 

The FWHM thickness of the HI disk was determined by fitting 
a Gaussian to each vertical slice across the galaxy major axis 
in each HI channel map, enabling the gas layer thickness to be 
measured independently on each side of the galactic centre. To 
facilitate vertical thickness measurements, each HI cube was ro- 
tated in the plane of the sky to align the galaxy major axis with 
the x axis, ie. cube axes (x,z, v), (see Sect. 6 of paper I for ro- 
tated HI cubes for each galaxy). The measured FWHM thickness 
is corrected for the effect of the beam by 



FWHM,,/// = ^FWHM* , c , - FWHMg 



(3) 



to obtain the actual gas thickness FWHM-,///, where the FWHM 
of the telescope beam is denoted by FWHM<j. For simplicity, 
hereafter we use FWHM- to refer to the beam-corrected HI gas 
thickness. 

The projected radius was deprojected to the galactocen- 
tric radius using the kinematics measured from the XV map 
in Sect. [2] to identify the galactocentric radius that contributes 
emission at the projected radius and line-of-sight velocity of the 
channel map being measured. This was facilitated by building a 
"face-on" HI cube that showed the HI distribution in the equa- 
torial plane (R', s) (where R' is the projected radius and s is the 
distance in the line-of-sight) of the emission in each observed 
line-of-sight velocity. For each galaxy a face-on HI cube was 
built using the radial kinematics determined in Sect. [2] 

To simulate the observational effects of the telescope, the 
face-on HI cube was convolved by the FWHM of the telescope 
beam along the x axis. In order to account for the intrinsic HI 
velocity dispersion and the finite spectral resolution, the spectral 
axis of the face-on HI cube was convolved by a ID point-spread- 
function with total velocity dispersion cr Vitot (R) equal to 



o- v ,tot(R) 



ct v ,hi(R) 2 + 



( FWHM„ ) 



2-y/21og2 , 



(4) 



where cr v j iI (R) is the radial HI velocity dispersion measured in 
Sect.|2]and FWHM y c /, is the spectral resolution of the telescope 
correlator. The telescope signal response function is a sine func- 
tion, which has a FWHM,,. e /, of * 1.2 times the width of the 
correlator channels. 
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Fig. 9. ESO074-G015 (IC5052). Top panel: Measured FWHM, thickness of the HI gas disks. Lower panels: Each of the kinematic 
functions of the HI reflected about the galactic centre to show both sides of the galaxy, and the average of both sides. Middle left: 
rotation curve; middle right: HI velocity dispersion; lower left: HI surface density; lower right: HI vertical thickness FWHM Z . 



The intensity-weighted velocity fields of the face-on HI 
cubes were calculated for each galaxy. These velocity fields pro- 
vide a representation of the galactocentric radii contributing flux 
to each major axis slice of the actual observed HI cubes. 

Each channel map of the face-on HI cube was masked by the 
maximum radius of HI emission on each side of the galaxy. Next 
the galactocentric radius at each projected radius of a channel 
map was derived by taking the intensity-weighted mean of the 
galactocentric radii of all points along that projected radius in the 
face-on channel map, where the intensity weights are the HI flux 
at those positions. This process was used to find the weighted 
mean galactocentric radius corresponding to the major axis po- 



sition of each vertical thickness measured from the actual ob- 
served HI cubes. 

The FWHM Z (/?) was determined from the average of all the 
FWHMj. fitted at that galactocentric radius. The number of chan- 
nel maps with emission from any one galactocentric radius usu- 
ally ranged from one to 20, as dictated by the galaxy kinematics 
and the resolution of the observations. Galaxies with flatter rota- 
tion curves enabled a greater number of measurements for each 
galactocentric radius. Higher spectral resolution relative to the 
HI velocity dispersion caused the HI emission at each galacto- 
centric radius to be spread over more channel maps. However, 
increasing the spectral resolution resulted in lower signal-to- 
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Fig. 10. ESO109-G021 (IC5249). Measured FWHM, thickness of the HI gas disk and the kinematic functions of the HI. For details 
see the caption of Fig. [9] 



noise and fewer fits, as fitting was only possible on vertical slices 
with at least 4 consecutive points with flux greater than two times 
the rms noise. 

To correct for the scatter in the vertical thickness measure- 
ments due to beam-smearing and low signal-to-noise, the whole 
HI cube fitting process was iterated in a similar manner as that 
used for the XV modelling algorithm. The known kinematics 
(from Sect. 13 and the newly measured FWHM 2 (7?) flaring fits 
formed the first model which was used to build a new Hl(x, z, v) 
HI cube. The effect of the telescope observation was simulated 
by convolving the HI cube spatially by the telescope synthesised 
beam and spectrally by total velocity dispersion 

Q~ v,tot 

(R) (de- 

fined in Eqn. ©). The model HI cube was then fitted, and the 



fits to the model HI cube compared to the fits of the actual ob- 
served HI cube. This cycle was iterated, with all kinematics ex- 
cept the flaring held fixed, until FWHM,(/?) fits to the model HI 
cube matched the FWHM»(/?) fits of the actual observed HI cube 
to within 10% of the FWHM of the telescope beam along the z 
axis or until the total difference in fits stopped decreasing. The 
best fits were usually obtained after 3-4 iterations. 

This iterative flare fitting method was tested under vari- 
ous conditions on synthetic galaxy observations. Simulations 
showed that high signal-to-noise observations could be corrected 
for beam smearing to within 5-10% of the FWHM beam size 
along z; however, lower signal-to-noise observations with large 
beams were somewhat limited with a FWHM- uncertainty of 
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Fig. 11. ESOl 15-G021. Measured FWHM, thickness of the HI gas disk and the kinematic functions of the HI. For details see the 
caption of Fig. [9j 



« 20% of the FWHM beam size along z. Although the flare 
fitting algorithm does accurately correct for the projection and 
beam-smearing of the HI distribution, the final FWHM- accu- 
racy is naturally still dependent on the adequate signal-to-noise 
of the HI data cube. 

4. Results for flaring profiles 

For each galaxy in our sample we present the best fit measure- 
ments of the FWHM- thickness of the HI disk as a function of 
radius, derived independently on each side of the galaxy. The 
measured thickness along each major axis slice of the HI chan- 
nel maps was undertaken using the kinematics measured by XV 



modelling in Sect. [2] The error bars shown on the flaring profiles 
in this section are the lcr errors of the Gaussian fits. 

We were able to measure the HI flaring to radii within 1-2 
kpc of the maximum radius at which the other HI kinematics 
were measured. This was a particularily impressive result con- 
sidering that the signal-to-noise of most vertical slices was only 
3-5. The average HI cube peak signal-to-noise was only 13.6, 
except for UGC7321, which had almost ~6 times the signal- 
to-noise of the other galaxies. Analysis of the vertical shape of 
the HI distribution in UGC7321, facilitated by the high sensitiv- 
ity observations, found unexpected azimuthal variations in the 
gas thickness and midplane displacement of this edge-on galaxy. 
These surprising results for UGC7321 are presented in Sect. 14.81 
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Fig. 12. ESO138-G014. Measured FWHM Z thickness of the HI gas disk and the kinematic functions of the HI. For details see the 
caption of Fig. [9j 



4.1. ESO074-G015(IC5052) 

The FWHM Z vertical thickness of ESO074-G015 expanded from 
FWHM Z = 0.4 kpc at a galactocentric radius of 1.5 kpc to 
FWHM Z = 2.5 kpc at « 7.0 kpc, with an additional thick- 
nesses of 3.3 and 4.0 kpc measured at the edge of the HI disk 
at radii 8.0-8.1 kpc (Fig. [9]). These outermost two measurements 
are consistent with the flaring shape of the disk, but the mea- 
surement has a large uncertainty of 0.5-0.6 kpc and were only 
measured in one channel map of the HI data cube. The flaring 
at 8 kpc comprises the warped eastward extension at a radius of 
4' in the HI column density map (see Fig. lb in paper I), and is 
detected as very diffuse edge emission in the outer channel maps 



on the receding side. The channel maps show the emission at 2 
x rms noise of 1.83 x 10 19 atoms cm -2 . 

The integrated flux of ESO074-G015 was 30% less than 
that measured from the HI Parkes All Sky Survey (HIPASS) 
line spectrum obtained from the online HIPASS Public Data 
Releas43- This suggested that there is additional extended HI 
emission that was not observed with the long ATCA array con- 
figurations. Indeed, a low resolution HI column density map 
available in the online Local Volume HI Survey (LVHIS0), 
showed large scale emission around the warped disk extend- 



1 http://www.atnf.csiro.au/research/multibeam/release 

2 http://www.atnf.csiro.au/research/LVHIS/ic5052/ATCAresults 
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ing to ~12 kpc. The smallest array configuration observations 
included in our analysis were low sensitivity snapshot observa- 
tions with the ATCA 375 array configuration. The LVHIS image 
was obtained with a 12 hour integration using the ATCA EW367 
array configuration. Additional low resolution data of ESO074- 
G015 may detect HI flaring to thicknesses greater than 4.0 kpc. 

Despite the gross asymmetry of the ESO074-G015 rotation 
curve, the flaring is remarkably symmetric (see Fig . O . This HI 
disk does thicken by -200 pc at radii inside 4 kpc particularly 
on the east side. This appears to be due to the increased HI ve- 
locity dispersion from ~6 to 8 km s~'at the same radii, and is 
more pronounced on the east side. The inner west side of the 
disk from 2-4 kpc shows intense star formation and radio con- 



tinuum emission, but the east side is optically more quiescent. 
The flaring profile is clearly well-resolved by the synthesised HI 
beam FWHM of 292 pc. However, the low peak signal-to-noise 
of 14 in the HI data cube restricts measurement of flaring in the 
faint outer HI disk. Despite being well-resolved spatially, the ap- 
parent mass asymmetry indicated by the rotation curve precludes 
this galaxy from dark halo flattening modelling. 

4.2. ESO109-G021 (IC5249) 

Like UGC7321, ESO109-G021 is renowned for being a "su- 
perthin" galaxy, exhibiting a scalelength-to-scaleheight ratio 
of 11. However, although ESO109-G021 has a similar 
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Fig. 14. ESO274-G001. Measured FWHM, thickness of the HI gas disk and the kinematic functions of the HI. For details see the 
caption of Fig. [9j 



maximum rotation speed to UGC7321, its radial scale is approx- 
imately 2.5 times larger, with a long scalelength of 7 kpc and a 
sharp stellar truncation at 17 kpc. We find that the HI disk flares 
from a FWHM Z of 800 pc at 2 kpc radius out to 2.8 kpc at 20 kpc 
(see Fig. [Toll. The thickness of the gas disk appears to be roughly 
constant from 2 to 5 kpc and begins to rise exponentially from a 
radius of 6-7 kpc. The steep flaring of the gas disk at 17-18 kpc 
is consistent with the edge of the stellar disk. 

The gas layer thick ness was also been measured by 
Ivan der Kruit et all (120011) from the total HI intensity map, as 
the signal-to-noise of the channel maps was too low as they had 
only the original 12 hour ATCA 6.0C observation by Carignan 
& Beaulieu. They found the FWHM, thickness rose from 1.3 



kpc near the centre to 4.9 kpc at 15 kpc. Our flaring measure- 
ments were obtained using an additional 30 hrs of ATCA time 
split over EW352, 750D and 6. OA array configurations, yielding 
a peak channel map signal-to-noise of 10.1 



The HI disk was possibly unresolved at radii less than 2 kpc, 
as the large synthesis beam of 1.18 kpc made it difficult to re- 
solve spatial features smaller than ~800 pc. An unresolved cen- 
tral gas thickness would explain the relatively thick central gas 
layer, comparable in thickness to the stellar disk. 
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Fig. 15. ESO435-G025 (IC2531). Measured FWHM, thickness of the HI gas disk and the kinematic functions of the HI. For details 
see the caption of Fig. [9] 



4.3. ESO115-G021 

ESO115-G021 shows strong flaring (see Fig. [Till from a cen- 
tral FWHM- thickness of 0.8 kpc to over 2 kpc at the outermost 
measured point, 6.1 kpc from the galactic centre. The HI emis- 
sion of the galaxy extends to a galactocentric radius of 8.2 kpc 
on the eastward side; however, this emission is very diffuse and 
too faint for HI thickness measurement. 

The peak signal-to-noise of the channel maps used to mea- 
sure the flaring was only 18.6 and the FWHM beam size was 194 
pc. Each radial fit was the average of all vertical thickness fits at 
major axis positions of a channel map that deprojected to this 
radius. Some thickness measurements were averaged from mea- 



surements to 5-7 channel maps; however, most were obtained 
from two to three channel maps. 

Despite the low average HI velocity dispersion of (cr v ) = 
5.5 ± 0.8 km s" 1 , the HI disk flares to a vertical HWFM Z above 
the plane of > 2 kpc near the HI edge. 

4.4. ESO138-G014 

ESO138-G014 was the only galaxy in our sample in which it was 
not possible to satisfactorily measure the HI gas layer thickness. 
The measured flaring thickness are shown in Fig. [T2l however, 
the derived FWHM- are ballpark estimates only due to miss- 
ing large scale structure that was not measured with the large 
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Fig. 16. UGC7321. Measured FWHM- thickness of the HI gas disk and the kinematic functions of the HI. For details see the caption 
of Fig. E 



ATCA array configurations used. ESO138-G014 and ES0274- 
G001 were the only two galaxies in our sample with extended 
HI disks (~10' diameter) for which we were unable to obtain 
compact array HI observations. Although they both had similar 
peak channel map signal-to-noise of ~15, the low spatial resolu- 
tion of ESO138-G014 prevented correction of the fitted flaring 
by iterative modelling. The irregular measured flaring profile, 
which displays a decreasing disk thickness at radii outside of 15 
kpc, is probably due to the missing large scale structure. 

Under a constant hydrostatic pressure, a fall in gas thickness 
of ~0.5-1.0 kpc could be explained by a decrease in HI velocity 
dispersion of 15-30% or 1.5-3.0 km s _1 . Although a 3 km s _1 
decrease in velocity dispersion is observed between 7 and 10 kpc 



radius, this decrease is significantly inside the onset radius of the 
decline in gas layer thickness. Due to the irregular gas flaring 
profile we decided not to model the dark halo of ESO138-G014. 

4.5. ESO146-G014 

The HI disk of ESO146-G014 flares dramatically from a 
FWHM, thickness of 0.4 kpc at radii of 3 kpc, to FWHM, = 1 .6 
kpc at 9.0 kpc from the galactic centre (see Fig. [T3l . The flaring 
is highly asymmetric around the galactic centre, probably due to 
the low peak HI cube signal-to-noise of 11.3 which limited the 
number of fits at each radius. Despite the noisy data, the mea- 
sured flaring profile still showed a steep exponential rise from 
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Fig. 17. Left - All the vertical profiles in a sample of galactocentric radial bins in UGC7321. The radius of the bin is given in the top 
right-hand corner of each plot. Each radial bin is 145.4 pc wide - the same size as the image pixels along the major and minor axes 
of UGC7321. The ordinate axis of each plot shows the HI column density in units of HI atoms cm" 2 in each vertical slice. Right 
- The same as on the left-hand side, but now we show average profiles, which were obtained after aligning the profile centres and 
vertically rescaling all profiles to a peak of 1. Overlaid in blue is the Gaussian fit to each average profile. 



4 kpc radius to the last measured point at 8 kpc radius. This 
steep inner flaring, occuring at 1.5 to 2.5 stellar scalelengths, is 
consistent with a low stellar surface density as inferred from the 
B-band face-on surface brightness of 23.8 mag arcsec" 2 and the 
low rotation speed. The low stellar surface brightness and large 
flaring make this galaxy an excellent candidate for measuring 
the shape of the dark halo. 

4.6. ESO274-G001 

As mentioned above, the HI imaging of ESO274-G001 is also 
missing some extended large scale structure due to a lack of 
HI observations using array configurations shorter than 1 .5 km. 
This missing structure was clearly apparent in the flaring mea- 
surement process as the measurements of the HI cube modeled 
with the initial measured flaring were much lower than the flar- 
ing measured from the observed HI channel maps. The best fit 
model displayed in Fig. [14] shows that the flaring rises roughly 
linearly from a FWHM. thickness of 200 pc at less than 1 kpc 
radius to 1.9 kpc at 5.5 kpc radius. 



4.7. ESO435-G025(IC2531) 

IC2531 displays a very large outer and remarkably symmetric 
flare (see Fig. [T3I ). rising from a gas thickness of 800 pc at 4.5 
kpc radius to 3-4 kpc at 30 kpc. The central gas thickness was 
initially partially unresolved given the FWHM,? beamwidth of 
1.3, and was decreased to the final thickness by iterating. The 
large gas thickness uncertainty of between 100 and 500 pc was 
due to the low peak channel map signal-to-noise of 10.9. 



4.8. UGC7321 

The measured flaring of UGC732 1 is our most accurate determi- 
nation of the flaring in our sample, rising from a FWHM. thick- 
ness of 350 pc at a radius of 0.5 kpc to roughly 4 kpc at 1 1 kpc 
from the galactic centre. The 16 hours of VLA time obtained 
by Matthews enabled the HI channel maps to be imaged with 
much better sensitivity than that acquired for the other galax- 
ies observed with the ATCA. As a result the error bars on the 
UGC7321 FWHM ; thickness measurements are only 15-50 pc, 
except in the outermost 2-3 kpc. The relatively flat rotation curve 
also allowed the flaring at most radii to be independently mea- 
sured in over 10-20 different channel maps. 

The flaring profile of UGC7321 exhibits the high degree of 
symmetry of the gas thickness about the galactic center as ex- 
pected for an undisturbed, isolated galaxy. At all radii the flar- 
ing profiles are symmetric about the galactic centre to within 
the bounds of the errors bars (lcr errors are AFWHM Z < 50 
pc within a radius of 9 kpc). As the HI velocity dispersion as 
a function of radius is also very symmetric, this shows that the 
total vertical force must also be symmetric. Given a symmetric 
luminous density distribution, this implies that the vertical force 
and density of the halo are symmetric about the galaxy centre. 

The most striking feature of the HI flaring profile of 
UGC7321 is the change in the shape of the flaring profile near 
the edge of the stellar disk. In the inner disk the HI thickness 
appears to increase linearly with radius, while outside the stellar 
disk it displays an exponential profile. In UGC7321, this "break" 
occurs at 7. 8 kpc, which i s betw een the stellar truncation radius 
of 6.7 kpc (Pohlenet al. , 120031) and the ma ximum stellar disk 
radius of 8.1 kpc dMattthews & Woodll2003l) . 
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Fig. 18. Measured HI FWHM Z thickness (upper panel) and HI 
midplane location (lower panel) of UGC7321 as a function of 
radius, where the measurements are colour coded by the line-of- 
sight velocity of each channel map. The most bluewards points 
were measured from the extreme approaching velocities, and 
those most redwards were measured from the extreme receding 
velocities. 
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The vertical thickness of the HI disk o f UGC7321 has also 
been studied by Mattthews & Wood (2003) who, as in the study 
of NGC2403 bv iFraternali et all i2Q0^K tried to discern the 
HI fla ring from the HI total intensity map. Matt thews & Woodl 
(2003) found that the best fit to the observed HI total intensity 
map was achieved using a thin HI disk that flares from a FWHM Z 
of 260 pc at the center to 1 .64 kpc at 11.8 kpc, and an addi- 
tional low density halo HI component of con stant FWHM 7 of 3.3 
kpc th at lagged behind the disk gas rotation. Matt thews & Woodl 
(2003) parameterised the thin HI disk using hf — h\ + hfor 5 ^ 4 . 
This produced a smooth flaring profile, but lacked the steep ex- 
ponential increase we measured in gas thickness at radii out- 
side of the stellar disk of UGC7321 and m ost of the galaxies in 
our sa mple (see Sect. 15.41 ). The addition by Mat tthews & Woodl 
(2003) of a low density second gas component with a large con- 
stant FWHM ; seems to simulate the thick gas layer at the outer 
radii. The filamentary nature of this structure suggests that it is 
HI ejecta, originating in the inner disk, rather than a pervasive 
diffuse thick HI layer. As shown below, our direct flaring mea- 
surements of UGC7321 show that the gas layer can be well rep- 
resented by a single gaussian gas component with a FWHM that 
varies smoothly with radius: modelling the gas layer thickness 
by two gas components appears to be unnecessary. 

The high signal-to-noise of these observations allowed in- 
spection of the intrinsic shape of the vertical gas thickness. The 
left-hand panel of Fig.[T7]displays the observed vertical profiles 
in a sample of galactocentric radial bins on both sides of the 
galactic centre. The profiles show the number of independent 
measurements in each bin and the flaring of the gas thickness on 



Fig. 19. Measured FWHM Z thickness (upper panel) and mid- 
plane displacement (lower panel) of the HI disk of UGC7321 
deprojected onto the X-Y plane, where the color scale on each 
plot shows the respective FWHM Z or midplane displacement 
measurements. Measurements in white are off the top end of 
the color scale. Each line of measurements forming a coloured 
track shows the measurements from a single HI channel map. 
The measurements are the same on each side of the Y — line, 
as due to the edge-on orientation of the gas disk, we can't de- 
termine if HI emission comes from the near or far side of the 
disk. 



both sides of the galaxy. In order to inspect the vertical shape 
more closely, we averaged the profiles in each bin after aligning 
the profile centres and vertically rescaling all profiles to a peak 
of 1 (see Fig. [17] right-hand panel). This shows that the profiles 
are indeed Gaussian in nearly all radii bins, with small high-z 
wings seen only in annuli at radii of ~4 and ~8 kpc. 

In Fig. [18] we show the direct FWHM Z (upper plot) and mid- 
plane position (lower plot) fits of the HI disk. The FWHM Z 
shows a clear increase in thickness at radii between 3.5 and 5.5 
kpc on both sides occuring on both sides of the galactic centre. 
The colors show that this occurs only in the extreme velocity 
channels, corresponding to HI in the line of nodes. The lower 
plot shows that the HI midplane also deviates at this radius and 
azimuth in the X-Y plane, dropping below ~200 pc below the 
average midplane position on the left-hand side, and rising ap- 
proximately 150 pc above the galactic plane on the right-hand 
side. 
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Fig. 20. The HI rotation curves for all galaxies in our sample overlaid on each other for comparison. 



In Fig.[19]we show the FWHM Z thickness (upper panel) and 
midplane position (lower panel) measurements deprojected onto 
the X-Y plane, where the color scale on each plot shows the 
respective FWHM Z or midplane displacement measurements. 
Each line of measurements forming a coloured "track" shows 
the measurements from a single HI channel map. The measure- 
ments are the same on each side of the Y=0 line as, due to the 
edge-on orientation of the gas disk, we can't determine if HI 
emission comes from the near or far side of the disk. Below we 
briefly describe the HI warp of the midplane before discussing 
the unusual azimuthal variation of both the thickness and mid- 
plane displacement mentioned above. 

The midplane displacement plot (lower panel) shows that 
disk plane warps by 300-500 pc in the plane of the sky on the 
left-hand side centred around Y=0, while on the right-hand side 
the warp amplitude is similar, peaking at a small angle of ~30° 
to the line-of-nodes. This implies that the measured thickness is 
uneffected by the warp on the left-hand side of the centre, while 
on the right-hand side the measured thickness of the direct loop 
1 fit may have been slightly overestimated due to the projection 
of the warp. These effects of warping are corrected by iterating 
the measurement as explained in Sect. 13.21 

Also visible is a pronounced deviation of the midplane and 
FWHM ; thickness of gas in the line of nodes at radii between 
2.5-5.0 kpc on the left-hand side of the galaxy and at radii from 
4-5.5 kpc on the right-hand side. This midplane displacement is 
seen most clearly in the lower plot of Fig. [18] while this devi- 
ation of the gas thickness from the mean flaring can be seen in 
both the upper plot of Fig. [18] and the upper panel of Fig. fl"9l 
One explanation for the increased broadening and tilting of the 
disk midplane along the Y=0 axis at this radius, is that the 
gas is responding to the local perturbation of the total poten- 
tial by the suspected bar. Peanut-shaped de viations from ellipti- 
cal isophotes in the projected R-band image dPohlen etal]l 2003) 
favour a bar oriented in the plane of the sky (ie. along the Y=0 



axis), rather closer to end-on which would align the bar axis with 
the line-of-sight. The measured enhanced local gas thickening 
along the Y=0 axis probably is the best observational evidence 
of a bar in an edge-on galaxyQ Other explanations for localised 
regions of increased thickness are infalling material from a com- 
panion or disruption by a companion, however UGC7321 is an 
isolated galaxy. The high symmetry of the flaring on both sides 
of the centre, and the high stellar disk major- to-minor axis ratio 
(a/b — 10.3 at fiR = 25.0 mag arcsec~ 2 , Matthews, Gallagher & 
van Driel 1999) render disruption by an undetected neighbour- 
ing galaxy rather implausible. 

Another noticable feature of Fig. |T9]is that the tracks denot- 
ing the location of the measurements in the X-Y plane display 
wiggles which systematically vary over a range of azimuthal an- 
gles. As our deprojection method is based on the velocity field, 
these wiggles could be due to deviations from circular rotation in 
the velocity field causing small errors in the deprojectioifl Kinks 
are typically seen in the velocity fi elds of less inclined galaxies 
with strong spiral arms (e.g. M81. lRoti . l 19741) suggesting that, 
if seen face-on, UGC7321 may exhibit weak spiral arms in the 
gas distribution. 



5. Discussion 

5.1. HI rotation curves 

In Fig. [20] the rotation curves of our galaxies are overlaid. IC 
2531 stands out as a massive, large galaxy compared to the rest 
of the sample. We have two very slow rotators with maximum 
rotation speeds of ~60 km/s in our sample, while the others form 
a group of galaxies rotating at ~100 km/s. The shape of the ro- 
tation curves are in good agreement with the rotation c urves of 
60 less inclined late-type dwarf galaxies measured by ISwatersI 



3 We thank Lia Athanassoula for pointing this out. 

4 Private communication with Agris Kalnajs. 
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(Il999h . Although low mass LSB galaxies exhibit a shallower in- 
ner rotation gradient, our high resolution observations clearly re- 
so lve differential r otati on, in agreement with oth er observations 
bv lSwatersi d 1999b and lde Blok & Bosmal d2002h . 

5.2. HI velocity dispersion 

The HI velocity dispersion measurements presented in this paper 
comprise the highest spatial resolution measurements of veloc- 
ity dispersion to date and the first measurements using data that 
span the radial extent of each galaxy by more than 10 HI syn- 
thesised beams. Fig. I2T1 shows that the HI velocity dispersion of 
our galaxies ranges from 4.5 km s _1 to 12.5 km s _1 , with most 
galaxies displaying a HI velocity dispersion of 6.5 to 7.5 km s _1 . 
Of the eight galaxies in our sample, all display radial structure 
in HI velocity dispersion, except ESO274-G001 which has an 
approximately constant HI velocity dispersion of 6.5 km s _1 . 

Fig. [22] shows the mean HI velocity dispersion as a func- 
tion of the maximum rotation velocity of the HI disk. There is 
a definite increase in the mean velocity dispersion with rotation 
speed from about 5 to 8 km s _1 for rotation velocities 70 to 120 
km s" 1 , but IC2531 has the same velocity dispersion of approxi- 
matel y 8 km s" 1 although i ts maximum rotation speed is 230 km 
s -1 . iDalcanton et all d2004l) found that edge-on galaxies with ro- 
tation velocity greater than 120 km s _1 have well-defined dust- 
lanes, while those with lower rotation do not. Dalcanton et al. 
suggested that well-defined dust lanes are related to lower gas 
turbulence and lower scale heights for the dust layer in the more 
rapidly rotating galaxies. Do our data support this suggestion? 

From Eqn. (1), the scaleheight h of a Gaussian dust layer is 
related to the ISM velocity dispersion cr by 



IT 



dK z 
' dz 



4nGp 



(5) 



where p a is the mid-plane total density which is dominated by 
baryons in normal spirals. The density p a = £/2/i* where 2 is 
the typical surface density of a disk and h, is the baryon scale 
height. It turns out that E and the stellar scaleheight are both ap- 
proximately proportional to the maximum circular velocity V c 
dGurovich etaTL 120091: iKregel et all l2002h . so the typical value 
of p is independent of V c and we expect the scale height of the 
dust layer to be directly proportional to the ISM velocity disper- 
sion cr. In our sample of galaxies, the galaxy with the largest V c 
(IC2531) shows a clean well-defined dustlane as expected from 
the Dalcanton et al. sample. However Fig. 22 shows no evidence 
for a decrease in cr as V c increases, so our data do not suppport 
the variable turbulence explanation for the change in dust lane 
morphology with V c . 

Two of the larger galaxes, IC2531 and ESO138-G014, dis- 
play HI velocity dispersions 2-3 km s _1 larger at small radii 
coinciding with the inner stellar disk. IC2531 is a large Sc 
galaxy with a strong bar, while the HI XV map of ES0138- 
G014 indicates a possible bar. The other barred galaxy in our 
sample, UGC7321, displays the next highest central HI veloc- 
ity dispersion. These three galaxies have a relatively large stan- 
dard deviation in the HI velocity dispersion as often observed 
in more face-on strongly barred galaxies like NGC1365, pre- 
sumably caused by gas shocks and increased star formation. 
Unfortunately, ESO138-G014 was not in our optical sample due 
to its low galactic latitude, so we were unable to confirm the 
existence of a bar by investigating the stellar isophotes. 

The velocity dispersion profile of UGC7321 exhibits three 
broad localised peaks with velocity dispersion 1-2.5 km s _1 
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Fig. 22. Mean HI velocity dispersion (cr v H1 (R)) against the max- 
imum rotation speed V max of each galaxy. The error bars dis- 
play the standard deviation of the HI velocity dispersion for each 
galaxy, showing that more massive galaxies tend to show larger 
radial variation of the HI velocity dispersion. 



higher than the intervening medium. These occur at the galactic 
centre, and at radii of 6 kpc and 10 kpc. Stellar observations of 
UGC7321 bv lPohlen et al.l (l2003h show a central peanut-shaped 
stellar over-density with a radius of 2.0 kpc, which could be 
caused by resonant off-plane thickening via bar buckling (e.g. 
Combes etal, 1990). The rise of HI velocity dispersion from 
6.5 km s _1 to 9.5 km s in this region suggests that processes 
associated with bar buckling could be providing the additional 
energy heating the gas. The second HI velocity dispersion peak, 
from 4 kpc to 7 kpc, corresponds to the radius of the bar shown 
in the "figure-8" HI sig nature and the stellar isophotes analysis 
by iPohlen et al.1 d2003l) . The third HI velocity dispersion is of 
lower amplitude and occurs outside the stellar disk. These re- 
sults suggest that at least part of the additional gas heating seen 
in the high inner HI velocity dispersions of barred galaxies in 
our sample could be due to star formation along radial flows in 
the bar, and bar buckling. 

Observations of the face-on galaxies NGC3938 

NGC628 
NGC1058 

19841 




van der Kruit & Shostak, 
Shostak & van der Kruit, 

(Ivan der Kruit & Shostak! Il984t ~lPetric & Rupert [2006), 
NGC6946 dBoulanger & Viallefondl 1 19921) and NGC5474 
(Row nd et all 11994) have shown that the HI velocity dispersion 
ranges from 4 to 16 km s _1 , and has generally been assumed 
to be 6-8 km s in the absence of signifi cant star formation. 
In NG C 69 46, detailed obs e rvatio ns by iKamphuis & Sancisil 
dl993l) and iBoomsma et al.1 d2008l) indicate the presence of 
widespread high velocity gas with significant kinetic energy, 
while the radial HI velocity dispersion profile shows a peak 
of more than 10 km s within R25, followed by a gradual 
decline. Large-scale turbulence probably plays a role in this. 
The cause of the HI velocity dispersions is complicated and 
not well-understood, but turbulence driven by star formation, 
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Fig. 21. The HI velocity dispersion profiles for all galaxies in our sample overlaid on each other for comparison. 



supernovae , shearing motions and magnetic ins tabilities all may 
play a role (iTamburro et an.l2009HAgertz et aUl2009l) . 

The typical velocity dispersion of 6.5 to 7.5 km s of galax- 
ies in our sample is in agreement with these observations, which 
is to be expected given the LSB nature of all the galaxies in our 
sample (except IC2531 and IC5052). On small scales it is as- 
sumed that increased HI velocity dispersion is due to heating 
of the ISM by stellar winds f rom hot OB stars, s upernovae and 
expan ding HII regions (e.g. Mc Kee & Ostriken Il977t ISpitzerL 
fl978h . Star formation is known to increase the turbulence of 
the ISM (lOort & Spitzerl 1 19551) . as the newly formed stars eject 
gas from the cloud within which they formed by radiation pres- 
sure and stellar winds. Although limited by a small sample size, 
our results suggest that increased gas velocity dispersion in disk 
galaxies is linked to the presence or prior formation of bars, 
which could cause the widespread star formation required to 
form strong dust lanes. However, it could also be caused by star 
formation in spiral arms. 

The most recen t observations of NGC1058 by 
Petric & Rupen (2006) are the only HI velocity dispersion 
study of a face-on galaxy with high spectral resolution (2.58 km 
s _1 ). They found a general radial trend of a declining HI velocity 
dispersion which was punctuated by many peaks and saddles 
where it ranged from 4 to 14 km s _1 . Two of the three regions 
with high HI velocity dispersion did not coincide with star 
formation tracers such as Hq, or radio continuum emission. For 
the third one it did, but that region was located at the galactic 
center. The high spatial frequency of HI velocity dispersion 
variation measured in our sample also suggests that the small 
scale velocity structure of the neutral ISM is likely to be quite 
complex. 

If the canonical HI velocity dispersion of 7 km s _1 found in 
the outer disks of galaxies were due solely to thermal motion, 
the gas temperature would need to be 7j « 6000 K. In the outer 
Galaxy where it is possible to measure the kinetic temperature 



and the velocity dispersion, Spitzei] (1 19781) found a thermal con- 
tribution to the HI velocity dispersion of only ~1 km s _1 , indi- 
cating that most of the broadening is non-thermal in origin. 

One possible explanation is that differential galactic rotation 
transforms rotational energy into turbulent motions, perhaps 
via shocks resulting fr o m disk instabilities dOuirogal 1983; 
iRohlfs & KreitschmannL [l987t ISellwood & BalbusL fl999b . 
Shocks in galactic bars could also generate gas heating. High 
resolution hydrodynamical simulations, including gaseous self- 
gravity, realistic radiative cooling, and rotation by IWada et all 
(120021) show that large scale rotation combined with gas self- 
gravity could be the source of the additional non-thermal energy 
broadening the gas velocity dispersion. Their models of turbu- 
lent rotating disks are consistent with the velocity dispersion 
of the extended HI disk in the dwarf galaxy NGC2915, which 
appears to have no prominent star formation. Their numerical 
simulations also show that the stellar bar in NGC 2915 enhances 
the HI velocity dispersion, which supports our finding of higher 
velocity dispersion in the inner disks of IC2531, ES0138-G14 
and UGC7321 

5.3. HI radial surface density 

Fig. [23] shows the measured radial HI surface density of each 
galaxy in our sample. UGC7321, ESO274-G001 and ESOH5- 
G02 1 all show a steep exponential decline at large radii, while 
IC5052 and ESO138-G014 appear to exhibit a roughly linear 
decline. With the exception of IC5052, most galaxies appear to 
show roughly smooth HI surface densities. 

5.4. Flaring of HI layers 

In Fig. [24] we display a plot of all the flaring profiles for each 
galaxy in our sample. The HI flaring profiles appear to obey 
a typical shape, with the flaring increasing linearly with radius 
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Fig. 23. The face-on surface density profiles for all galaxies in our sample overlaid on each other for comparison. 



where the stellar disk dominates the local gravitational potential, 
and steepening to an exponential profile in the outer disk where 
the potential is probably dominated by the halo. This break in 
the change of the flaring profile is seen most clearly in UC7321. 
The onset of exponential curvature appears to occur just outside 
the stellar truncation radius in the more massive galaxies with 
maximum rotation speed V max > 100 km s _1 . In smaller, lower 
surface brightness disk galaxies with V max < 100 km s _1 , the 
break possibly occurs at smaller radii due to dark matter domi- 
nating the inner gravitational potential. 

All of the galaxies in our sample display this flaring pro- 
file, except ESO138-G014 which is considered unreliable due to 
missing large-scale HI structure (see Sect. 14.41 ). Unfortunately, 
we were unable to detect this radial change in flaring profile 
shape in ESO274-G001, the galaxy with the highest spatial res- 
olution in our sample, as we were unable to measure the gas 
thickness at radii outside the stellar disk, due to a lack of short- 
spacing observations. Excluding ESO138-G014, all the more 
massive galaxies (V ma x > 100 km s~') display a shallower flaring 
profile over the inner disk. With the exclusion of ESO146-G014, 
the three remaining small disk galaxies {V„ mx < 100 km s" 1 ) all 
display steep linear inner flaring profiles. 

The flaring profile of ESO146-G014 shows remarkable sim- 
ilarities to the flaring profile of UGC7321. Apart from a similar 
central face-on B-band stellar surface brightness (respectively 
23.4 and 23.8 mag arcsec~ 2 ), these two galaxies appear quite 
different — ESO146-G014 has a much lower rotation speed and 
HI velocity dispersion. 

The three well-measured large galaxies (IC253 1 , IC5249 and 
UGC7321) and the unusual ESO146-G014, each display one 
common profile with a small FWHM Z . A common flaring pro- 
file in the inner disk suggests a common ratio of the HI veloc- 
ity dispersion to the stellar density, assuming the stellar density 
dominates at low radii. The ratio of mean HI velocity dispersion 
to maximum rotation speed is very similar in IC5249, UGC7321 



and ESO146-G014, the value for IC2531 is roughly half that of 
the others. 

The radial gradient of the flaring is more varied in the 
smaller galaxies. The coalignment of the inner flaring slopes of 
UGC7321, IC5249 and IC2531 is probably a coincidence, as 
a comparison of IC2531's flaring profile to that of the Galaxy 
(lOlling & Merrifieldll2000T) shows that the gas disk of IC2531 is 
approximately twice as thick as the Galactic HI layer. 

In Fig. [25] we show the flaring profiles with the horizontal 
axis scaled with (half) the diameter of the HI distribution, taken 
from paper I. This does make the profiles more comparable in 
shape, although still there is quite some variation of the thick- 
ness of the HI in kpc from galaxy to galaxy. We have done the 
same by scaling to (half) the optical radius, taken from NEE0- 
Not surprisingly, the result is similar since the optical and HI di- 
ameters do in general correlate among galaxies, although the rel- 
ative positions of the lines with respect to each ot her is different . 
For c omparison, the flaring profile of NGC 891 (Ivan der KruitL 
1981) would be on the scale of Fig.|25]run from about 0.3 kpc at 
0.2 scaled galactocentric radius to 1.9 kpc at 1.0. This is thinner 
than in our systems, but that is not surprising since our sample 
was selected to be made up of galaxies with low contributions of 
the disks to the gravitational potential and the thickness of the HI 
layer is to a larger extent determined by that of the dark matter 
halo, which would be less flattened than the disk. 

In this work we have measured the vertical thickness of the 
HI disk by fitting a single Gaussian gas component. We have 
not attempted to fit the thickness by a two Gaussian compo- 
nents as others have done, because the single Gaussian model 
was a good fit. Earlier higher resolution VLA observations of 
NGC891 yielded an exponentially flaring thin HI layer. More 
recent studies of its gas thickness, using exceptionally sensitive 



5 This NASA/IPAC Extragalactic Database is operated by the Jet 
Propulsion Laboratory, California Institute of Technology, under con- 
tract with the National Aeronautics and Space Administration. 
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lower resolution WSRT observations with w 350 hrs integra- 
tion time, recover extended filamentary high latitude HI with 
unusual velocity structure comprising 20% of the total HI mass 
dOosterloo et all 12001 fFraternali et"alll2004 . By applying this 
dual gas component mo delling proces s to t he initial 144 hrs of 
WSRT data on NGC89 l. lSwaters etail (1 19971) found that the tra- 
ditional warped and flared gas layer did not simulate the high lat- 
itude inner disk gas. Indeed the best models of NGC891 tested 
in this fashion used two disks each of constant thickness, with 
FWHM Z;(/ll „(fi) ~ 0.5 kpc and FWHM z , ftet (7?) ~ 6 kpc, for the 
thin and thick disk respectively, where the thick disk rotation 
lags that of the planar gas, and also displays some radial infall 
motion. 

Similar two-component HI modelling has been under- 
taken on seve ra l less -inclined larg e galaxies: NGC240 3 
dFraternali et all 120021). N GC4559 (tBarbieri et all l2005h 
NGC253 dBoomsmal 120071) and NGC6946 dFraternali et all 



2004; Boomsmal l2007h . Each of these galaxies is actively star- 
forming with significant fractions of their total HI mass (10- 
30%) exhibiting anomalous vel ocities forbidden by circula r ro- 
tation. As discussed in Sect.|H Mattth ews & Wood! ((2003) also 
applied this model to UGC7321, a low surface brightness galaxy 
in our sample. In each of these galaxies, modelling has found a 
thick gas layer combined with a thin constant thickness or flaring 
thin gas layer. With the exception of UGC7321, the high stellar 
density, star formation rate and stellar disk area of these galax- 
ies are all conducive to ejecting a greater HI mass fraction to 
high latitudes, whether by discrete star formation or large scale 
galactic fountains. Despite the massive amount of HI in the ha- 
los of these galaxies, this inner high latitude HI is probably more 
akin to Galactic high velocity clouds than a continuous broad HI 
component. 

As discussed in paper I, a small amount of high latitude fila- 
mentary HI is observed in the larger galaxies in our sample, with 



V max > 100 km s _1 , at projected radii corresponding to the inner 
disk. However this gas makes up only a small fraction (<sc 5%) 
of the total gas mass and may not be in hydrostatic equilibrium. 
For this reason we have not attempted to measure its mass dis- 
tribution. 



5.5. Open issues 

As mentioned above, the HI velocity dispersion measured in 
face-on galaxies were all undertaken with low spatial resolution 
of 1 kpc, much larger than the scale of most processes in the ISM 
which are likely to cause local variation in the HI velocity dis- 
persion. Deep high spatial and spectral HI observations of such 
systems are lacking. 

In theory, given high spatial and spectral resolution, it is pos- 
sible to determine the intrinsic shape of the HI line spectrum, and 
whether the velocity structure is best rep resented by one or more 
Gaussians, or a non-Gaussian function. iPetric & Rupenl {2006) 
found that the HI line profile was best fit by a two-Gaussian 
shape where the ratio of the integrated flux of the thick and thin 
Gaussians was constant at all radii with good signal-to-noise. 
The only edge-on galaxy with bright HI that is close enough 
to image at the requ ired re s olution is the small Scd galaxy, 
NGC4244. studied bv lolirn"gldl996al) . 

It is also important to test the isothermality of the HI distri- 
bution by fitting the HI velocity dispersion to XV maps formed 
from gas at different heights above the galactic plane. This is im- 
portant; without further information about the z-dependence of 
the velocity dispersion, we would need to assume z-isothermality 
when using the equations of hydrostatic equilibrium to derive 
the vertical potential of the dark matter halo. Although the spa- 
tial resolution of ESO274-G001 is sufficient to independently 
measure the gas at several heights above the plane, additional HI 
observations in compact array configurations would be needed. 
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Our HI observations in paper I show significant optical depth 
(assuming HI spin temperatures of 200-400 K) in isolated loca- 
tions of several galaxies in our sample. It would be interesting 
to generalise the XV modelling to also fit the z direction, as this 
would allow the HI opacity to be accurately measured across the 
galaxy. Such a generalised model would measure the gas den- 
sity p(R, z). By integrating over the line of sight, the total HI col- 
umn density and HI opacity consistent with the kinematic model 
could be calculated and compared to the observed HI column 
density to check the calculated HI opacity. This requires higher 
spatial resolution data than we have available in this study. 

6. Summary and conclusions 

Using our new XV modelling method detailed in paper II we 
have successfully fit the rotation curves of all eight galaxies in 
our sample; and the HI velocity dispersion and surface density in 
all the galaxies except for the inner disk of IC253 1 . The rotation 
curve was the most reliably fit achieving an accuracy of roughly 
2-4 km s _1 . The surface density was also very reliably measured. 
The HI velocity dispersion was more sensitive to local variations 
and noise in the XV map, however despite these challenges the 
uncertainty in the derived HI velocity dispersion is only about 1— 
1.5 km s _1 , which is compa rable to the accuracy obtained from 
NGC4244 by 011in3 dl996ah . The prevalence of bars in the larger 
galaxies and the increase in the mean velocity dispersion with 
maximum rotation speed suggests that bar buckling may be a 
common cause of widespread gas heating in disk galaxies. 

The measured rotation curves show that all galaxies in our 
sample display differential rotation, even the smallest galaxies 
IC5052 and ESO115-G021 which have peak rotation speeds of 
only 60 km s _1 . The increased velocity dispersion by 1-2.5 km 
s _1 observed in galaxies with bars suggests that shocks caused 
by disk instabilities in bars can cause gas heating. 

Iterative fitting of the HI layer thickness enabled us to reli- 
ably measure the HI flaring in 7 of the 8 galaxies in our sample, 
despite FWHMg beamwidths of up to 1.3 kpc and a low peak 
channel map signal-to-noise of around 13 for most of the galax- 
ies in our sample. The only galaxy we were unable to reliably fit 
was ESO138G-014 which was missing large scale HI structure 
due to an absence of short baseline observations, and also suf- 
fered low peak signal-to-noise (15) and a large telescope beam 
(0.965 kpc). For galaxy observations with only two of these three 
handicaps it was possible to iteratively recover a flaring model 
that fit the observed HI channel map distribution. UGC7321 is 
the most accurately measured flaring profile in our sample with 
errors of only 15-50 pc over most of the radial extent. 

The flaring profiles show an impressive degree of the sym- 
metry between flaring measurements on either side of the galaxy 
— this is most clearly seen in UGC7321. Axial symmetry of 
the gas layer thickness is expected as in the absence of local 
variations of the stellar surface density or gas velocity disper- 
sion, the vertical flaring should reflect the smooth symmetric 
gravitational potential of the dark matter halo. The exception- 
ally smooth flaring profile of UGC7321 shows that, on the 0.7 
kpc scale of telescope beam, the stellar surface density and dark 
matter halo gravitational potential are also smooth. 

The most remarkable feature of HI flaring profiles is a com- 
mon shape displaying a linear increase at low radii and steep- 
ening to an exponential increase in the outer disk. The steep- 
ness of the inner flaring should be dominated by the radial sur- 
face density, while at the stellar disk truncation radius the flaring 
becomes exponential and is dominated by the dark halo shape. 
This typical flaring shape is found in all of the galaxies except 



ESO138-G014 which was incompletely imaged, and ES0274- 
G001 which appears to have a more linear flaring profile. With 
the exception of ESO146-G014, all of the less massive galaxies 
(V max < 100 km s _1 ) display steep inner flaring profiles. The four 
larger galaxies (V max > 100 km s" 1 ), excluding ESO138-G014, 
display much shallower inner flaring profiles. 

The high sensitivity of the UGC7321 observations allowed 
us to discern variations from azimuthal symmetry of both the 
HI midplane position and the HI thickness. We show that at 
radii from approximately 2.5 to 6 kpc along the Y=0 axis, the 
FWHM, gas thickness varies by up to 500 pc from the az- 
imuthally averaged flaring in those radial annuli. At this same 
locations along the Y=0 axis the gas midplane tilts, dropping 
~200 pc below the average disk midplane on the left-hand side, 
and rising ~200 pc above the midplane on the right-hand side. 
These local deviations are considered to be best observational 
evidence for a bar in an edge-on disk galaxy. 

The next paper in this series will be devoted to an analysis 
of the data on UGC7321 as a first attempt at deriving constraints 
on the flattening of a dark halo. 
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